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The size of two-body weakly bound objects : short versus long range potentials 
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The variation of the size of two-body objects is investigated, as the separation energy approaches 
zero, with both long range potentials and short range potentials having a repulsive core. It is shown 
that long range potentials can also give rise to very extended systems. The asymptotic laws derived 
for states with angular momentum £=1,2 differ from the ones obtained with short range potentials. 
The sensitivity of the asymptotic laws on the shape and length of short range potentials defined by 
two and three parameters is studied. These ideas as well as the transition from the short to the long 
range regime for the i = 0 case are illustrated using the Kratzer potential. 


The study of the properties of weakly bound systems has 
found a renewed interest after the discovery of halo nuclei 
in nuclear physics [jlj -ftl as well as loosely bound dimers 
in molecular physics [l|.(ij . These systems have very large 
mean square radii and small separation energies. They 
can be treated as two (or three)-body systems interact¬ 
ing through a potential. In fact, the separation energy of 
the one or two nucleons forming the halo is so small that 
their degrees of freedom can be separated from those of 
the nucleons constituents of the core. In the case of dif¬ 
fuse diatomic molecules, the situation is even better, the 
separation energy being several orders of magnitude less 
than the ionization potential of each constitutent. 

In this Letter, we want to discuss in very general terms 
the size of two-body weakly bound objects. Short and 
long range repulsive forces have been studied up to now 
in the context of nuclear and molecular physics. In par¬ 
ticular, the variation of the size of these quantum systems 
as the energy approches zero (asymptotic laws) has been 
actively investigated [^||7| . 

There is one more class of potentials which can give rise 
to very extended systems in the zero energy limit, namely 
long-range potentials for which this limit is attained when 
the deepness of the potential goes to zero. It is the pur¬ 
pose of this paper to discuss these long-range potentials 
and show that they predict a variation of the size of the 
system, in the zero energy limit, different from short- 
range potentials. We will compare the two behaviours 
and discuss the transition between the short and long 
range regime. 

Concerning short range potentials, we will discuss 
short range potentials presenting a repulsive core. The 
presence of a repulsive core in the interacting potential 
should actually be included to modelize the Pauli exclu¬ 
sion principle. We will show how the presence of the core 
modifies the asymptotic laws found for short range po¬ 
tentials, discussed up to now, which are defined by only 
two parameters, namely the deepness of the potential and 
the range. 

The Letter is structured as follows. We will first show 
that asymptotic laws for states of any angular momentum 
l can be derived starting from the Schrodinger equation, 


without making any specific reference to the potential. 
Then we will discuss how we can derive these laws by us¬ 
ing the asymptotic behaviour of the wave functions both 
for short and long range potentials. Finally we will use 
the Kratzer potential to illustrate the energy dependence 
of the mean square radius for long range potentials, the 
influence of the core in the case of short range poten¬ 
tials (obtained from the Kratzer potential by cutting it 
at some distance) and the transition from the short to the 
long range regime. We will use the results obtained for a 
square well for a comparison with short range potentials 
defined by only two parameters. 

In this work, we will consider only spherical symmetric 
potentials. The spin degrees of freedom are ignored. In 
this case, the radial Schrodinger equation (Ti = 2m = 1) 
is given by 

[--x~2 ~ “7T + ^ X ^ + 'W r )]'M r ) = E e ^i(r). (1) 

or z r or r z 


where A scales the deepness of the potential w(r). Since 
we consider only the lowest states of each angular mo¬ 
mentum (the wavefunction has no node), they are simply 
labelled by I. 

The first very general result obtained from eq.(!jf), with¬ 
out any specific reference to the potential, is a kind of 
Heisenberg relation ]sj : 


< r 2 > t > 


(2 1 + 3) 2 
4 <T> e ’ 


( 2 ) 


relating the rms radius to the average kinetic energy. For 
confining potentials, we have Ee >< T >i, the equality 
being reached as A —> 0. Moreover, for power-law poten¬ 
tials (w(r) = r“, a > —2) (or superpositions of them) 
the virial theorem 


< T > e = ^ J pe{r)r[^-]d 3 r . (3) 


ensures that < T >goc< V >t and thus < T >^oc E£. 
Consequently, for this large class of potentials, the in¬ 
equality (2) readily tells us that the rms radius is diverg¬ 
ing with 1/E{ for all £ as E( —> 0, which means A —> 0. 
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No such a general prediction can be made in the case 
of short range potentials, because in this case the zero 
energy limit is obtained when the coupling A tends to a 
critical finite value A c |J. 

In the case £ = 0, another prediction is given by 
the Bertlnrann-Martin inequality [fil| which concerns the 
ground state radius and yields 


< r 2 > 0 < 


3 

E\ — Eq 


( 4 ) 


where Eq and E\ are the energies of the ground and first 
excited states, respectively. Noting that E\—E 0 > 0, we 
can write Ei = ipEo and get 


< r 2 > 0 < 


3 

(<P - l)^o 


ip — 1 > 0 . 


( 5 ) 


Consequently, for £ = 0 the asymptotic behaviour of the 
rnrs radius as the energy tends to zero is given by this 
expression for all, short or long range, potentials. 

In order to specify the value of ip and obtain predic¬ 
tions for £ ^ 0, it is necessary to introduce assumptions 
concerning the wavefunctions. 

Let’s first discuss short range potentials, characterized 
by a finite range Rq. Since the potential becomes negligi¬ 
ble beyond Rq, the asymptotic form of the wavefunctions 
is related to the spherical Hankel functions 


p-IJ-r 

Mr) « — . (6) 

As p —> 0 the energy tends to zero. Therefore, as long as 
the contribution from the inner part of the wavefunction 
is negligible, the rms radius is given by 

/* oo /»oo 

< r 2 >t& / e~ 2tir r 2 ~ 2t dr / / e~ 2tir r~ 2e dr . (7) 

J Rq J Ro 


With this approximation, we get the relations derived by 
Riisager, Jensen and Moller |ll| : 


< r >o" 


Co 2 .. _ 

^ !<r >1? 


Ci Ro 


VWl 


; < r 2 > 2 ~ c 2 Rq , 


( 8 ) 


Intuitively, the above assumptions imply the constants 
eg, to be independent of the potential. However this state¬ 
ment is not exact and in general ce = ci{w) is expected. 
For £ > 3 the rms radius behaves in the same way as 
for £ = 2, namely it tends to a constant as the energy 
approaches zero. 

Coming back to the inequality (j|), in the case the £ = 1 
state is in the continuum, Ei is set to zero (or tp = 0). 
(Note that this procedure is not valid for confining poten¬ 
tials, for which E( > 0.) Using the asymptotic proper¬ 
ties of the s-state wavefunction leads to an absolute lower 
bound (l2) 


< r 2 > 0 > 


2|£o| ' 


The equality is reached as —Eq —> 0, independently on 
the potential. As a consequence, in eq.(|j), cq = 1/2 is a 
firm prediction, independently of the shape of the short 
range potential. 

To get a practical insight on how much C£ depends on 
w, the asymptotic laws (||) have also been studied numer¬ 
ically for potentials depending only on two parameters, 
the deepness and a typical length Rq. In particular, the 
square well and the gaussian potentials have been used 
]lT| . Their typical lengths are the square well radius and 
the range of the gaussian, respectively. In ]llj] , it is shown 
that the coefficients ce seem to be not very sensitive to 
the particular shape of the potential. 

In general, < r 2 > 1 = f(E,Ro,w), that is the mean 
square radii depend on the energy of the state, the typ¬ 
ical length R 0 and the shape of the potential. Eq.(||) 
explicitly give the dependence of < r 2 >£ on Rq for 
short range potentials defined by two parameters. This 
can be seen directly from the Schrodinger equation, by 
means of the change of variable x = t/Rq. As a result, 
< x 2 >£= f'(e,w) where e = EiRq. As soon as the 
potential contains more than two parameters, the scal¬ 
ing cannot be exact, although it may constitute a good 
approximation under some circonstances. As a conse¬ 
quence of this gwasf-scaling, the < r 2 >£ may depend on 
the length of the potential in a way similar to (|j). Such 
a situation may occur in particular when the potential 
has a short range repulsive component. 

Let us now discuss the asymptotic laws for long range 
potentials, using the behaviour of the wave functions. 
Contrary to short range potentials, the wavefunction is 
in this case confined inside the potential. Use can be 
made of the spherical Bessel function of the first kind 
ji(kr), with k —> 0, cutting the integrals at the first zero. 
This procedure yields right away 


< r 2 >*« 


ct{w) 

\Ee\ 


( 10 ) 


The derivation is very crude, and merely confirms the 
above statement of eqs.(||) and ini- 

Relations (|), (||) and @ show that both short and 
long range potentials may predict very diffuse systems 
in the zero energy limit. The asymptotic laws obtained 
present similarities as well as differences. In fact, in the 
case of long range potentials, the rms diverges as 1/Et 
not only for £ = 0, as it has been found for short range 
potentials, but for any £. Besides, the constants C£ are 
sensitive to the potential w. 

A question one may ask about systems with halos is 
how to define them. In the case of short range potentials, 
a natural definition arises as it has been extensively dis¬ 
cussed in the literature 0. A reference scale is needed 
to compare it to the mean square radius of the system in 
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order to quantify the size of the halo. This reference scale 
can be taken equal to the range of the short range poten¬ 
tial which in some cases can be related to a physical scale 
of the system, for example the size of the core in a halo 
nucleus and is identified with the outer classical turning 
point. An ideal halo is then defined as a system for which 
the single halo particle has a very large probability to be 
inside the classically forbidden region. Its properties are 
therefore determined by the tail of the wavefunction and 
are almost independent of the potential. 

For the s-state, the above discussion refers to E 0 — > 0. 
However, it is interesting to note that for the shell-delta 
potential, w(r) = S(r — Rq) which possesses a single 
bound s-level, the particle is always in the classically for¬ 
bidden region for any finite Eq. 

In the case of long range potentials, the definition of a 
halo seems less straightforward. If the potential presents 
a repulsif core, a natural reference scale is given by the 
size of this core which can be related again to some phys¬ 
ical size as, for example, that of the core in a halo nucleus 
or the sum to the two atomic radii, in the case of a diffuse 
molecule. So, the mean square radius of the system can 
be compared to the range of the potential core in order 
to quantify the extension of the halo. Contrary to the 
case of short range potentials, the wavefunction of the 
single halo particle will always be confined in the clas¬ 
sically allowed region and the probability of finding the 
particle in the potential will always be one. Finally, in 
the case of short range potentials, an ideal halo is only 
possible if the halo particle is in an s-state. In fact, for 
p- a d-states the centrifugal barrier plays an important 
role, confining the wavefunction. This is not the case 
for long range potentials, where the barrier plays no role 
and therefore wavefunctions of any angular momentum 
can be extended. 

Let us now discuss the asymptotic laws for long and 
short range potentials as well as the transition between 
the two in a particular case, namely the Kratzer potential 

@ : 

w(r) = —2A(~ - | 1 ) (11) 

where a is just a scaling parameter which gives the dis¬ 
tance at which the potential changes from attractive to 
repulsive. This potential has a long range and a repul¬ 
sive core. The short range case is obtained by putting 
the potential to zero at a variable distance R cu t '■ 
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Wcut{r) = —2A( — - )0(i? cut - r) (12) 

r Zr* 

This cut Kratzer potential has the peculiarity of being 
defined by three parameters, namely the size of the core, 
the deepness of the potential A and the range R cu t- 

By solving the Schrodinger equation (jlj) with ©, w e 
obtain the mean square radius for any t : 


{r 2 )z = 


m\ 


3 + 5 


7 2 + (^) 2 + 2( 7 2 + (^) 2 ) 


(13) 


where q 2 


(a 2 A). This gives, when Ei —> 0, 


(r 2 )o 



<r 2 )i 


7.5 

~E{ 



(14) 


As expected, ( r 2 )c have the same dependence on the en¬ 
ergy as predicted by (J|) and ©• For the sake of com¬ 
parison, note that the values of ci are the same in the case 
of a pure Coulomb force. For the harmonic oscillator, Co 
and ci take the same values whereas C 2 = 10.5. 

We have solved |I]) with the cut Kratzer potential (0) 
for different values of R cu t in order to see explicitly how 
the asymptotic laws vary when we use a three instead 
of a two parameter potential. Particular attention has 
been devoted to the asymptotic region, working down to 
energies of the order of 10 -10 . Note that such an achieve¬ 
ment is only possible for potentials admitting analytical 
solutions for which the continuity conditions can be eas¬ 
ily solved with high accuracy, whereas solving numeri¬ 
cally the Schrodinger equation is very uncertain at low 
energies. The calculations we present are performed for 
o=l. The features we emphasize in this work are not 
qualitatively affected by this parameter. From the nu¬ 
merical results, we get the following approximate asymp¬ 
totic laws : 


Here C\ and C 2 slightly depend on R cut , whereas Co = 1/2 
as expected ( 

Relations (15|) show that, for any £, the dependence of 
the mean square radii on the energy of the state is not 
modified even in presence of a core. This is in agree¬ 
ment with what the arguments based on the asymptotic 
behaviour of the wave functions would suggest (|j). 

Concerning the dependence of the mean square radii 
on R cut , we see from (|l5|) that the cut Kratzer potential 
present a gitasi-scaling, the dependence on the range of 
the potential R cut being very similar to that of (||). 

In order to show the dependence of the q on R cu t and 
the shape of the potential, these coefficients are shown, 
as a fonction of the energy, in figure 1 (up t = 0, middle 
l = 1, bottom l = 2), for R cut = 3 (line), 6 (short-dashed 
line) and 12 (long-dashed line) . For comparison we show 
the ci obtained for a short range potential defined by two 
parameters (dotted line). We take the square well as an 
example. 

We see that, apart from the £ = 0 case, as the en¬ 
ergy tends to zero, the asymptotic values of ci depend 
indeed on the potential. This dependence increases with 
£, as one should expect if the centrifugal barrier has a 
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confining effect and the wavefunctions become more sen¬ 
sitive to the particular shape of the potential. We also 
see that for 2 = 1,2 the spreading among the three cut 
Kratzer potential results is only of a few percents (legiti¬ 
mating the gwasi-scaling), much less than the differences 
with the square well potential results. This qualitatively 
underlines the differences in the inner part of the wave- 
functions. 

It is interesting to ask if specific physical systems are in 
the low-energy limit where the asymptotic laws directly 
give the size of the system once the separation energy is 
known. We have looked at the case of halo nuclei and 
molecular dimers, by using experimental values of the 
separation energy and radius to calculate the correspond¬ 
ing co- For 11 Be the experimental separation energy is 
.5 MeV and the rms radius is about 7 fm, which leads to 
Co « 1.1. The 4 He — 4 He dimer has a separation energy 
of 10 _1 /reH and a radius of 70 A Q ; it yields Co = .52. 
The calculated 3 He — 6 Li [|j, with a separation energy of 
10 _2 /reK and a radius of 210 A, reaches Co = .5. These 
estimates show that the loosely bound dimers are situ¬ 
ated in the asymptotic region, whereas the halo nuclei 
merely remain at its threshold, in a range sensitive to 
the potential. A similar conclusion has been drawn by 
studying the Bertlmann-Martin inequality |p~4|| . 

Outside the asymptotic region, it becomes even diffi¬ 
cult to decide between a short and a long range potential, 
if the only observables are the separation energy and the 
rms radius. The key point in this respect is obviously 
given by the spectrum. In both dimers and halo nuclei, 
the number of states is rather limited, and thus a descrip¬ 
tion by a short range potential is quite natural. 

Let us now come to the transition between the short 
and long range regime. The example of the Kratzer po¬ 
tential provides us with a way of studying this transition 
by looking at evolution of ct as a function of the range 
Rent ■ A priori, the transition depends on the energy, but 
numerical results show that this dependence is simple in 
our case. In figure 2 we present co against \J\Eo\R cut . 
The half-point between the short range (co = 1/2) and 
long range regime (co = 3) lies at around ^/\Eo\R cu t = 2. 

For higher angular momenta, the problem is more com¬ 
plicated because going from the short © to the long 
range © regime modifies both ce and the energy de¬ 
pendence of the asymptotic laws. 

Note that the scaling and quasi- scaling present in 
eqs.(j|) and © cannot be used to study the transition 
between short and long range regime. For example, for 
t = 0, the transition does not occur if we take the square 
well potential at any finite radius Rq. This is due to the 
fact that according to eq.(U|J) co = 1/2 independently on 
Rq. Nevertheless, it is clear that for a fixed enegy, the re¬ 
quired potential deepness A is decreasing as Rq increases, 
so that the limit Rq —s► oo implies A —* 0. It means that 
the square well potential only reaches the long range limit 
for infinite Rq. 


Let us finally remark that up to now the weakly bound 
objects found experimentally are intepreted in the two- 
body description as s-states. It remains a challenging 
problem to observe halo states that would be interpreted, 
in this context, as higher angular momentum states. To 
the extent that a two-body description yields the basic 
description, they should occur as excited levels. In the 
nuclear case, this situation is hardly possible, because 
the many-body degrees of freedom will strongly affect 
the first s-state, and destroy the simple two-body picture. 
From this point of view, the dimers are better candidates. 

In conclusion, we have discussed the size of two-body 
weakly bound objects using both long range potentials 
and short range potentials having a repulsive core. We 
have shown that not only short range but also long range 
potentials give rise, in the zero energy limit, to very ex¬ 
tended systems. The asymptotic laws obtained in this 
case show a different energy dependence from those of 
short range potentials already for p- and d-states. The 
dependence of the asymptotic laws on the particular 
shape and typical length of short range potentials, de¬ 
fined by both two and three parameters, has also been 
studied. All these ideas as well as the transition of the 
asymptotic laws from the short to the long range regime 
for £ = 0 are illustrated in the particular case of the 
Kratzer potential. 
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FIG. 1. cut Kratzer potential : the values of ci are plotted 
against Et for i = 0 (up), 1 (middle) and 2 (bottom). The line 
corresponds to R cu t = 3; the short-dashed and long-dashed 
lines correspond to R cu t = 6 and 12, respectively. For compar¬ 
ison, the dotted line shows the results obtained with a short 
range potential without the repulsive core. As an example, 
the square well is taken. As far as Co is concerned, the arrows 
indicate the range of the lowest energies corresponding to halo 
nuclei (arrow on the right) and dimers (arrow on the left). 



FIG. 2. Evolution of Co from the short to the long range 
regime for the cut Kratzer potential. For comparison, the 
dashed line shows the value of 1/2 valid when Eo —> 0 for any 
short range potential (Jjj). 
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